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ANB  THE  LAW  BE  THE  WAKE 

by 

D.  M.  Nelson 


Underwater  Ordnance  Department 


ABSTRAC T.  A  method  for  calculating  the  incompressible 
turbulent  boundary  layer  based  on  the  ’’law  of  the  wall'*  and 
Coles*  ’’law  of  the  wake"  is  presented.  The  method  is  ap¬ 
plicable  to  two-dimensional  bodies  and  to  bodies  of  revolu¬ 
tion  in  axisymmetric  flow  when  the  boundary- lay er  thickness 
is  not  necessarily  small  compared  to  the  body  radius.  The 
carrying  out  of  this  method  involves  a  simultaneous  solution 
of  the  momentum  integral  equation  and  the  energy  integral 
equation*  assuming  that  the  mean  velocity  profiles  are  given 
by  a  universal  two-parameter  representation  as  suggested 
'  y  Coles . 

Since  these  two-parameter,  mean  velocity  profiles  inher¬ 
ently  contain  sufficient  information  to  determine  the  shear 
stress  at  the  wall  and  the  rate  of  energy  dissipation  in  the 
boundary  layer,  and  accurately  represent  velocity  profiles 
to  the  point  of  separation,  it  is  anticipated  that  there  will 
result  an  improved  technique  for  estimating  the  growth  of 
a  turbulent  boundary  layer  which  may  be  valid  to  the  point 
of  separation. 


s. 


NAVAL  ORDNANCE  TEST  STATION 
China  Lake,  Californla/^"^  ^ 


November  1964 


r\ 


oo 


K 

V  • 


o-  f>0 


u.  s. 


NAVAU  OltONANCK 


TKST  STATION 


AN  ACTIVITY  OF  THE  BUREAU  OF  NAVAL  WEAPONS 


J.  I.  Hardy,  Caft.,  USN 

Commander 


Wm.  B.  McLean,  Ph.CX 

Tmchnicat  Director 


FOREWORD 

This  report  presents  the  theoretical  development  of  a 
method  for  calculating  the  incompressible  turbulent  boundary 
layer  based  on  the  'Taw  of  the  wall"  and  the  "law  of  the  wake." 
This  development  was  carried  out  to  provide  a  more  rigorous 
solution  of  the  boundary- laye r  equations  for  turbulent  flow,  in 
order  to  obtain  an  improved  technique  for  estimating  the  inflow 
velocity  profiles  to  propellers  and  the  skin- friction  drag  of 
bodies  of  revolution. 

The  work  was  done  at  the  U.S.  Naval  Ordnance  Test  Station 
between  November  1962  and  January  1964  under  Bureau  of 
Naval  Weapons  Task  Assignment  RU TO- 3E- 000/ 21  6- 1 /F008 - 
06-02  Problem  Assignment  231.  H.  E.  Eggers  was  the  Bureau 
of  Naval  Weapons  Project  Engineer, 

The  considered  opinions  of  the  Propulsion  Division  are 
presented  in  this  report. 
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NOMENCLATURE 
Profile  parameter  (Eq.  13) 

An  empirical  constant  in  the  "law  of  the  wall"  (Eq.4) 
Unit  vector 

Shape  parameter  (Eq  .  7) 

Scale  factor 
Thickness  ratio,  ^r/^r 
Thickness  ratio,  6j*/0r 

Curvature  of  body  surface  in  plane  of  mean  flow 

Mixing  length 

A  length  defined  by  Eq.  27 

Numerical  constant  equal  to  2.75 

Numerical  constant  equal  to  2.50 

Numerical  constants  (see  Appendix) 

Order  of  magnitude 
Static  pressure  of  the  fluid 
Mean  static  pressure 
Fluctuation  static  pressure 
Radial  coordinate  (Fig.  3) 

Body  radius 
Reynolds  number 

Components  of  fluid  velocity  in  direction  of  increasing 
X,  y»  V  *  respectively 

Mean  values  of  fluid  velocity  components 
Fluctuation  values  of  fluid  velocity  components 
Fluid  velocity  at  edge  of  boundary 
Undisturbed  free- stream  velocity 

Velocity  of  fluid 

"Law  of  the  wake"  function 
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X,  y,  y  Orthogonal  curvilinear  coordinates  (Fig.  3) 

X  A  length  used  to  define  a  Reynolds  number 

P  Angle  between  tangent  to  body  surface  and  axis  of  symmetry 
6  Boundary- layer  thickness 

6^  See  discussion,  page  11 

6^  Displacement  thickness 

Pseudo-displacement  thickness  for  l/r^  /  0;  displacement 
thickness  for  l/r^  =  0 

Pseudo- energy  thickness  for  l/r^  /  0;  energy  thickness 
for  l/r^  =  0 

A  Turbulent  thickness  (Eq.  6) 

6  Momentum  thickness 

Pseudo-momentum  thickness  for  l/r^  /  0;  momentum 
thickness  for  l/r^  =  0 

«  Proportionality  between  mixing  length  and  distance  from 
the  wall 

V  Kinematic  viscosity 
n  Profile  parameter  (Eq.  5) 
p  Fluid  density 

T  Shear  stress 

Shear  stress  at  the  wall 
Laminar  shear  stress 
Turbulent  shear  stress 
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INTRODUCTION 

Coles'  boundary- layer  investigations  (Ref.  1)  showed  that  the  two- 
dimensionalt  incompressible,  turbulent  boundary-layer,  mean  velocity 
profiles  in  the  whole  range  from  transition  to  separation  closely  ap¬ 
proximate  a  two-parameter  family.  This  fact  suggests  the  possibility 
of  using  these  profiles  in  place  of  the  more  common  power-law  pro¬ 
files^  in  computing  the  development  of  the  turbulent  boundary  layer. 
Such  a  method  might  be  expected  to  yield  satisfactory  results  up  to  the 
point  of  separation,  whereas  power- law  profile  methods  tend  to  break 
down  near  separation  partly  because  at  this  point  the  mean  velocity 
profiles  cannot  be  accurately  approximated  by  a  power  law. 

Since  one  of  the  parameters  governing  the  two-parameter  family  is 
V#  =  n/ Tq/p ,  the  shear  stress  at  the  wall  is  specified  once  the  mean 
velocity  profile  is  specified.  Hence,  a  method  that  uses  these  two- 
parameter  velocity  profiles  would  not  have  to  make  use  of  an  approxi¬ 
mate  expression,  such  as  a  flat-plate  formula,  for  the  shear  at  the  wall. 
These  profiles  also  allow  the  rate  at  which  energy  is  dissipated  in  the 
boundary  layer  to  be  calculated.  In  general,  their  use  permits  a  more 
rigorous  approach  to  the  boundary-layer  solution.  Furthermore,  an 
investigation  based  on  available  experimental  data  (detailed  below)  in¬ 
dicates  that  the  mean  velocity  profiles  on  bodies  of  revolution  in  axi- 
symmetric  flow  also  fall  reasonably  well  w*ithin  this  two-parameter 
family  so  that  such  a  method  may  be  applied  almost  equally  well  to  a 
body  in  axisymmetric  flow. 

Using  Coles'  two-parameter  velocity  profile  family  as  a  starting 
point,  the  method  presented  here  for  determining  the  incompressible 
turbulent  boundary  layer  on  two-dimensional  bodies  and  bodies  of  rev¬ 
olution  in  axisymmetric  flow  was  undertaken.  Since  the  equations 
governing  the  axisymmetric  flow  about  a  body  of  revolution  reduce  to 
those  for  the  flow  over  a  two-dimensional  body  as  ro  -►  oo  (r^  =  body 
radius),  only  the  axisymmetric  development  is  given.  By  setting 
l/r^  =  0  at  any  point  in  the  axisymmetric  development,  the  equivalent 
two-dimensional  result  is  obtained. 


^  Power-law  profiles  are  of  the  form  u/U  =  (y/6)^. 
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HISTORICAL  REVIEW 


In  1925,  Prandtl  (Ref,  2)  made  an  important  contribution  to  the 
theory  of  turbulent  flow  in  the  development  of  his  mixing-length  hy¬ 
pothesis.  For  the  case  of  near-parallel  two-dimensional  flow,  this 
hypothesis  results  in  a  relationship  between  the  shearing  stress,  the 
mixing  length,  and  the  gradient  of  the  velocity  normal  to  the  direction 
of  the  flow. 


T 


(1) 


Applying  this  result  to  the  flow  near  a  wall,  Prandtl  assumed  that 
the  mixing  length  varied  proportionally  with  the  distance  from  the  wall 
and  that  the  shearing  stress  was  independent  of  the  distance  from  the 
wall. 

i  =  icy 

(2: 

T  =  Tq  =  pv^ 

The  first  of  these  assumptions  follows  from  the  fact  that  the  mixing 
length  must  go  to  zero  at  the  wall.  The  second  assumption  is  reason¬ 
able  if  the  pressure  gradient  in  the  direction  of  flow  is  small,  since 
for  near-parallel  two-dimensional  flow. 


1  ap  ax 

- ar  — 

p  ax  dy 


Substituting  Eq.  2  into  Eq.  1  and  integrating  from  u  =  0  to  u  =  u,  there 
is  obtained 


u  1 

-  =  —  (Iny  -  In  y^)  (3) 

ic 

where  y^  is  the  unknown  distance  from  the  wall,^  where  Ti  =  0.  Using  a 
dimensional  argument,  y^  must  be  proportional  to  v/v^.  Hence  Eq.  3 
takes  the  following  form: 

u  1  yv^ 

-  =  —  In  -  +  C  (4) 

V^  K 

This  is  the  logarithmic  mean  velocity  distribution  which  is  the  ‘Taw  of 
the  wall'*  for  two-dimensional,  incompressible  flow. 


^  Since  Eq.  3  is  not  applicable  to  the  region  very  near  the  wall  where 
the  laminar  stresses  predominate  (the  laminar  sublayer),  it  does  not 
satisfy  the  boundary  condition  u  =  0  where  y  =  0. 
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The  validity  of  Eq.  4  for  flow  past  a  wall  was  substantiated  by  the 
experiments  in  1940  of  Schultz-Grunow  (Ref.  3).  In  these  experiments, 
the  mean  velocity  profiles  and  skin  friction  were  measured  for  turbulent 
flow  over  a  flat  plate  with  no  pressure  gradient.  The  data  from  the 
inner  part  of  the  boundary  layer  was  found  to  be  in  good  agreement  with 
Eq.  4. 

An  additional  result  of  these  experiments  was  that  the  deviation  of 
the  velocity  profiles  from  Eq.  4  in  the  outer  part  of  the  boundary  layer 
could  be  expressed  quite  accurately  as  a  function  of  y/6  only.  Thus 
the  mean  velocity  profiles  for  turbulent  flow  over  a  flat  plate  with  no 
pressure  gradient  could  be  expressed  as 

u  1  yv^  /yv 

K  V  \  6  / 

The  fact  that  the  Schultz-Grunow  experiments  were  carried  out  with 
no  pressure  gradient  made  them  consistent  with  the  assumption  used  in 
developing  Eq.  4  that  the  Br/dy  is  very  small.  The  question  therefore 
arose  whether  Eq.  4  woula  also  be  valid  for  boundary  layers  developing 
in  a  strong  pressure  gradient  where  the  dr/By  would  no  longer  be  small. 
This  question  was  answered  by  the  experiments  in  1949  of  Ludwieg  and 
Tillman  (Ref.  4),  in  which  the  mean  velocity  proliles  and  skin  friction 
were  measured  for  turbulent  flow  over  a  flat  plate  in  both  positive  and 
negative  pressure  gradients.  Once  again  the  data  from  the  inner  part 
of  the  boundary  layer  (excluding  the  laminar  sublayer  and  transition 
zone)  were  found  to  be  in  good  agreement  with  Eq.  4.  However,  the 
deviation  of  the  velocity  profiles  from  Eq.  4  in  the  outer  part  of  the 
boundary  layer  was  found  to  be  no  longer  only  a  function  of  y/6  but  also 
a  function  of  the  distance  along  the  plate,  x.  Hence,  the  mean  velocity 
profiles  for  turbulent  flow  over  a  flat  plate  with  a  pressure  gradient 
could  be  expressed  as 


u  1  yv^  jy  \ 

—  =  -  In  -  +  C  +  F  -  ,  X  j 

K  V  \6  / 

In  1956,  Coles  (Ref.  1)  investigated  the  form  of  the  function  describ¬ 
ing  the  deviation  from  Eq,  4  of  the  velocity  profiles  in  the  outer  portion 
of  the  boundary  layer.  He  examined  essentially  all  the  experimental 
data  pertinent  to  the  growth  of  two-dimensional  turbulent  boundary 
layers  in  a  pressure  gradient,  and  determined  that,  to  a  very  close 
approximation,  this  function  of  y/6  and  x  could  be  expressed  as  the 
product  of  a  function  only  of  x  and  a  function  only  of  y/6.  Thus  Coles 
expressed  the  mean  velocity  profiles  for  two-dimensional  turbulent 
boundary  layers  in  the  following  form: 

u  1  yv^  n(x)  /y\ 

-  =  In  -  +  C  +  -  Wj-  j  (5) 

tc  V  K  \  s  I 


3 


NAVWEPS  REPORT  8510 


¥ 

•  Coles  refers  to  the  function  W(y/6)  as  the  'Taw  of  the  wake"  since  it 

represents  the  portion  of  the  flow  in  the  boundary  layer  that  is  charac* 
teristic  of  wake  flow,  i.  e,  ,  essentially  unaffected  by  Tq.  The  function 
W(y/6)  is  given  in  Table  1.  The  profile  parameter,  11,  together  with 
determines  the  shape  of  the  mean  velocity  profiles. 

TABLE  1.  Coles*  Wake 
Function 


y/6 

W 

y/6 

W 

0.00 

0.000 

0.55 

1.152 

0.05 

0.004 

0.60 

1.307 

0.10 

0.029 

0.65 

1.458 

0.15 

0.084 

0.70 

1.600 

0.20 

0.168 

0.75 

1.729 

0.25 

0.272 

0.80 

1.840 

0.30 

0.396 

0.85 

1.926 

0.35 

0.535 

0.90 

1.980 

0.40 

0.685 

0.95 

1.999 

0.45 

0.50 

0.838 

0.994 

2.000 

THE  TURBULENT  THICKNESS  AND  SHAPE  PARAMETER 

The  fitting  of  analytical  expressions  for  velocity  profiles  to  experi¬ 
mental  data  is  often  done  by  requiring  equal  momentum  thicknesses,  0  , 
and  equal  displacement  thicknesses,  5*.  However,  for  an  expression 
of  the  form  of  Eq.  5  it  ir  much  more  convenient  not  to  work  with  0  and 
6*  directly  but  to  make  use  of  the  turbulent  thickness,  A,  and  the  shape 
parameter,  G,  utilized  by  Clauser  (Ref.  5). 


It  can  easily  be  shown  that  two  velocity  profiles  having  equal  values  of 
A,  G.  and  v^  also  have  equal  values  of  6  and  6^. 
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DETERMINATION  OF  v*  FROM 
MEAN  VELOCITY  PROFILE 

Because  of  the  apparent  universal  applicability  of  Eq.  4  to  the  inner 
part  (excluding  the  laminar  sublayer  and  transition  zone)  of  the  turbulent 
boundary  layer,  it  may  be  used  to  make  a  reasonable  estimate  of  the 
value  of  v^ /U  if  only  the  mean  velocity  profile  is  known,  Eq.  4  may  be 
written  as  follows; 


U  1  1 

—  - - In  y  +  -  —  In  -  +  -  C 

U  U  «  U  K  V  U 


Taking  the  derivative  with  respect  to  In  y, 


V*  d(u/U) 

U  d(ln  y) 


(8) 


By  plotting  u/U  versus  In  y  and  determining  the  slope  of  the  linear  (or 
near-linear)  inner  portion  of  the  profile,  the  value  of  v^ /U  is  obtained. 


APPLICABILITY  OF  COLES*  RESULT  TO  THE 
BOUNDARY  LAYER  ON  A  BODY  OF 
REVOLUTION  IN  AXISYMMETRIC  FLOW 

The  boundary- layer  equations  for  the  axisymmetric  flow  over  a  body 
of  revolution  reduce  to  the  equations  for  flow  over  a  two-dimensional 
body  as  6/rQ  0.  Hence,  the  velocity  profiles  over  the  forward  part 
of  a  body  of  revolution  in  axisymmetric  flow  where  6/rQ  <<  1  should  be 

almost  identical  to  those  found  on  a  two-dimensional  body.  However, 
over  the  after  part  of  the  body  where  h/r^  is  no  longer  small,  the  de¬ 
viation  of  the  velocity  profiles  in  the  outer  portion  of  the  boundary- 
layer  from  the  law  of  the  wall,  Eq.  4,  might  be  expected  to  depend  on 
8/ro  in  such  a  fashion  that  Coles'  result,  Eq.  5,  for  two-dimensional 
flow  might  not  apply.  To  test  whether  Eq.  5  could  adequately  represent 
the  mean  velocity  profiles  for  the  axisymmetric  case,  the  measured 
mean  velocity  profiles  (Ref.  6,  7,  8)  on  the  after  part  of  several  bodies 
of  revolution  were  compared  with  Eq.  5, 

The  comparison  was  carried  out  in  the  following  manner.  Each 
profile  was  plotted  as  u/U  versus  In  y  and  an  estimate  of  v»/U  was  de¬ 
termined  from  Eq.  8,  This  value  of  v^ /U  and  the  curve^  of  7i/U  versus 


Due  to  the  lack  of  reliable  data  very  near  the  wall,  the  measured 
profiles  were  extended  analytically  to  the  wall  using  Eq.  4  and  the  equa¬ 
tion  for  the  flow  in  the  laminar  sublayer,  u/v*  =  yv^/v. 
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DETERMINATION  OF  v*  FROM 
MEAN  VELOCITY  PROFILE 

Because  of  the  apparent  universal  applicability  of  Eq.  4  to  the  inner 
part  (excluding  the  lanninar  sublayer  and  transition  zone)  of  the  turbulent 
boundary  layer,  it  may  be  used  to  make  a  reasonable  estimate  of  the 
value  of  Vy. /U  if  only  the  mean  velocity  profile  is  known.  Eq.  4  may  be 
written  as  follows: 


u  v^  1  1  v^^  V* 

—  - - In  y  +  -  —  In  -  +  -  C 

U  U  «  U  K  V  U 


Taking  the  derivative  with  respect  to  In  y, 


d(u/U) 

-  =  K  - 

U  d(ln  y) 


(8) 


By  plotting  u/U  versus  In  y  and  determining  the  slope  of  the  linear  (or 
near-linear)  inner  portion  of  the  profile,  the  value  of  /U  is  obtained. 


APPLICABILITY  OF  COLES'  RESULT  TO  THE 
BOUNDARY  LAYER  ON  A  BODY  OF 
REVOLUTION  IN  AXISYMMETRIC  FLOW 


The  boundary- layer  equations  for  the  axisymmetric  flow  over  a  body 
of  revolution  reduce  to  the  equations  for  flow  over  a  two-dimensional 
body  as  S/r^  0,  Hence,  the  velocity  profiles  over  the  forward  part 
of  a  body  of  revolution  in  axisymmetric  flow  where  h/r^  <<  1  should  be 
almost  identical  to  those  found  on  a  two-dimensional  body.  However, 
over  the  after  part  of  the  body  where  S/r^  is  no  longer  small,  the  de¬ 
viation  of  the  velocity  profiles  in  the  outer  portion  of  the  boundary- 
layer  from  the  law  of  the  wall,  Eq.4,  might  be  expected  to  depend  on 
6/ro  in  such  a  fashion  that  Coles*  result,  Eq.  5,  for  two-dimensional 
flow  might  not  apply.  To  test  whether  Eq.  5  could  adequately  represent 
the  mean  velocity  profiles  for  the  axisymmetric  case,  the  measured 
mean  velocity  profiles  (Ref.  6,  7,  8)  on  the  after  part  of  several  bodies  ^ 
of  revolution  were  compared  with  Eq.5. 

The  comparison  was  carried  out  in  the  following  manner.  Each  i; 

profile  was  plotted  as  u/U  versus  In  y  and  an  estimate  of  v^^./U  was  de-  r; 
termined  from  Eq.  8.  This  value  of  v^ /U  and  the  curve^  of  Ti/U  versus  S 


Due  to  the  lack  of  reliable  data  very  near  the  wall,  the  measured 
profiles  were  extended  analytically  to  the  wall  using  Eq.  4  and  the  equa¬ 
tion  for  the  flow  in  the  laminar  sublayer,  u/v*  - 
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y  were  used  in  Eq.  6  and  7  to  determine  values  of  A  and  G  by  numerical 
integration.  Applying  the  boundary  cond'  .ion  u  =  U  when  y  =  6  to  Eq.  5, 
there  is  obtained 

u  1  6v  n 

-  =  -  In  -  -I-  C  +  —  W(  1)  (9) 

Subtracting  Eq.  5  from  Eq.  9  yields, 


U  -  u 


1  y 

—  In  — 
«c  6 


+ 


n 

« 


W(l)  -  w 


)l 


(10) 


Substituting  Eq,  10  into  Eq,  6  and  7  and  carrying  out  the  integration,^  the 
terms  involving  W(y/6)  being  integrated  numerically  using  Table  1, 
there  results 


K  A 

6  = - 

1  +  n 

i.sizn^  +  (3.200  -  «G)n  +  2  -  kg  =  o 


Equations  11  permit  values  of  6  and  FI  to  be  determined  from  the  ex¬ 
perimentally  obtained  values  of  A  and  G.  Putting  these  values  of  6  and 
n  into  Eq.  5  rewritten  in  the  following  form. 


where 


a  direct  ccmparison  on  the  basis  of  u/U  versus  y  can  be  made  between 
the  experimental  profiles  and  Eq.  5.  The  comparisons  are  shown  in 
Fig.  la  -  d.  The  agreement  is  sufficiently  good  that  application  of 
Coles*  result  for  two-dimensional  flow,  Eq.  5,  to  a  body  of  revolution 
in  axisymmetric  flow  appears  reasonably  sound. 


Since  Eq.  5  or  Eq.  10  does  not  hold  in  the  region  very  close  to  the 
wall  (the  laminar  sublayer  and  the  transition  zone),  the  terms  involving 
ln(y/6)  give  rise  to  improper  integrals.  These  improper  integrals 
are,  however,  convergent  and  the  deviation  of  Eq.  5  or  Eq.  10  from  the 
true  flow  near  the  wall  can  be  shown  to  have  an  insignificant  effect  on 
the  values  of  parameters  such  as  A  and  G. 
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AN  ANALYTICAL  APPROXIMATION  TO  EQUATION  5 


Coles'  expression  for  the  mean  velocity  profiles,  Eq.  5,  has  two 
drawbacks.  First,  it  is  not  an  analytical  expression  since  W(y/6)  is 
given  as  a  table  of  values  so  that  manipulations  involving  integration 
must  be  carried  out  in  part  by  numerical  integration.  Second,  it  does 
not  satisfy  the  boundary  condition  at  the  edge  of  the  boundary  layer 
that  3u/8y  -  0  when  y  “  6.  Although  neither  of  these  drawbacks  seri¬ 
ously  impedes  its  usefulness,  it  would  be  desirable  to  have  a  simple 
analytical  expression  amenable  to  straightforward  integration  that 
would  closely  approximate  Eq,  5  and  also  satisfy  the  boundary  condition 
d\x/dy  -  0  when  y  =  6.  Such  an  v*:cpression  was  found  in  the  form 


u 

1 

A 

B 

fy 

V-+ 

=  —  In 

K 

+  C  +  — 1 

K 

b) 

b) 

Applying  the  boundary  condition  3u/0y  =  0  when  y  "  6,  Eq  .  \Z  becomes 

u  1  A  /  y  v"  1/1  n  \IV\^ 

_  =  _  In  -  +  C  +  --  -__  +  _A  -  (13) 

V)<.  K  V  tc  \m  m  I  \  6  I 


The  profile  parameter  A  in  Eq.  13  plays  the  same  role  as  the  profile 
parameter  O  in  Eq.  5,  By  choosing  n  =  2.50  and  m  =  2,75  in  Eq.  13,  it 
very  closely  approximates  Eq.  5.  A  comparison  between  Eq.  13  and 
Eq.  5  showing  this  fact  can  be  made  as  follows. 

If  Eq.  13  and  5  are  both  fit  to  the  same  velocity  profile,  the  boundary- 
layer  thickness,  6,  defined  by  Eq.  13  will  not  be  exactly  equal  to  that 
defined  by  Eq.5,  To  avoid  confusion,  from  this  point  on,  a  subscript  c 
will  be  put  on  the  boundary-layer  thickness  defined  by  Coles*  Eq.  5  to 
differentiate  it  from  the  boundary-layer  thickness  defined  by  Eq.  13, 

Applying  the  boundary  conditon  u  =  U  when  y  -  5  to  Eq.  13  and  sub¬ 
tracting  Eq.  13  from  this  result  yields 


Substituting  Eq.  10  into  Eq.  6  and  7  and  also  Eq.  14  into  Eq.  6  and  7, 
carrying  out  the  integrations,^  there  is  obtained 


(  14) 


and 


The  same  comments  made  in  footnote  4  concerning  Eq.  5  and  10 
apply  to  Eq.  13  and  14. 
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1  +n 

=  6e  - 


Z  +  3.^0Cn  +  1.52211^ 
G  = - 

«{ 1  +  ri) 


0.7333  +  0.04762A 
A  =  6  - 

K 

1.406  +  0. 1  162  A  +  0. 00  3297  A^ 

G  =  - 

*(0.7333  +  0.04762A) 


n  =  2.50 
m  =  2.75 


Specifying  that  both  profiles  have  the  same  values  of  A  and  G,  the  equa¬ 
tions  above  allow  the  determination  of  A  *  f(n)  and  b / ti.  -  ~  f(n). 

With  these  relationships  known,  a  direct  comparison  of  Kq.  5  and  13 
can  be  made.  In  order  to  make  this  comparison  independent  of  the 
value  of  V:j^ ,  Eq.  5  and  13  are  rewritten  in  the  following  form, 

u  1  Av4f  1  y  n  /  y  \ 

- In  -  -  c  =  -ln—  +  —  Wj -  (15) 

v^ic  V  K  A  tc  \6^} 

“u  1  A  v^,.  1  y  A/yV^  ^  \/y 

-  -  —  In  -  -  C  -  —  In  —  +  — I—  j  -  —  I —  +  —  Aj  j  —  j  (  16) 

K  V  K  A  «\6/  K  \m  m  f  \  6  I 

The  comparison  between  Eq.  15  and  16  is  shown  in  Fig.  2a  -  g  for  values 
of  n  ranging  from  0.1  to  5.0,  a  range  that  covers  the  values  most  com¬ 
monly  encounterea  in  practice.  The  agreement  is  excellent,  making  it 
reasonable  to  use  Eq.  13  and  14  in  place  of  Eq.  5  and  10  in  all  subsequent 
analyses.  Although  it  is  understood  that  n  =  2.50  and  m  =  2.75,  the 
symbols  n  and  m  are  retained  in  the  following  work  to  simplify  the 
writing  of  expressions. 


COORDINATE  SYSTEM 

The  most  convenient  coordinate  system  for  analyzing  the  axisym- 
metric  flow  over  a  body  of  revolution  is  the  orthogonal  curvilinear 
system  shown  in  Fig.  3.  Here,  x  is  the  distance  measured  from  the 
nose  along  the  surface  of  the  body  in  a  meridian  plane;  y  is  the  distance 
from  the  surface  of  the  body  measured  along  an  outward  normal  to  the 
body  surface;  y  is  the  angular  position  of  points  around  the  body;  K 
denotes  the  curvature  of  the  body  surface  in  the  meridian  plane;  and 
p  denotes  the  angle  between  a  tangent  to  the  body  surface  and  the  axis 
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FIG.  3.  Curvilinear  Coordinate  Systenn  Used  in  Analysis. 


of  symmetry.  The  scale  factors^  for  this  curvilinear  system  corre¬ 
sponding  to  the  coordinates  x,  y,  and  y  are 

h^  =  1  +  Ky 

hy  =1  (17) 

hy  =  r  =  r^  +  y  cos  p 


VECTOR  OPERATIONS 

In  an  orthogonal  curvilinear  coordinate  system,  the  operations  in¬ 
volving  the  del,  V,  operator  take  on  a  form  diff<*rent  from  that  t'n- 
countered  in  a  rectangular  system.^  For  an  orthogonal  curvilinear 
system  with  coordinates  xji  X3  and  scale  factors  h|,  h^.  h3,  the 

follov.  ing  expressions  are  valid.’ 


A  discussion  of  curvilinear  coordinates  may  be  found  in  Ref.  9. 
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B  =  div  B  = 


hih2h3 


V  X  B  =  curl  B  = 


1  d4> 

1 

8(t> 

1 

84> 

—  - e 

1  +  — 

02  + - 

ij  8x1 

9X2 

1^3 

8x3 

d 

8 

.(h2h3Bi)  +  - 

-(h3hiB2)  + 

- (hih2B3) 

8x2 

8x3 

hjei 

h3e3 

1 

9 

9 

9 

hih2h3 

9X1 

9x2 

8X3 

hiBi 

h2B2 

h3B3 

(18) 


In  the  above,  is  any  scalar  and  B  is  any  vector  having  components 
Bl*  B3  in  the  xj,  x2f  x3  system,  and  ei,  ez,  and  e3  are  the  unit 

vectors  in  the  directions  of  increasing  x^,  X2i  and  X3,  respectively. 
It  is  of  coarse  understood  that  the  xj,  X2f  X3  system  is  right-handed. 


THE  NAVIER-STOKES  AND  CONTINUITY  EQUATIONS 


The  equations  governing  the  flow  of  an  incompressible,  viscous 
Newtonian  fluid  are  the  Navier-Stoke s  and  continuity  equations.  For 
this  incompressible  case,  and  neglecting  gravitational  forces,  these 
equations  may  be  written  (Ref.  10)  as  follows. 


8V  ^  ^1^—1  ^ 

- V  X  (V  X  V)  +  -  V  (V  •  V)  =  Vp  -  vv  X  (V  X  V) 

dt  2  p 


V  .  V  =  0 


(19) 


In  these  equations  V  is  the  vector  velocity  of  the  fluid  and  p  is  the 
static  pressure.  The  components  of  V  in  the  directions  of  increasing 
X,  y,  and  \  will  be  denoted  as  u,  v,  and  w,  respectively. 

Expanding  Eq,  19  using  Eq.  17,  and  using  Eq.  18  with 


X]  =  X 

h]  =  h 

X2  =  y 

h2  =  h 

X3  =  Y 

h3  =  h 

there  is  obtained 
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(Flow  along  body) 
3u  u 


du  du 

—  +  V -  + 


w  du 
K  + - 


ar 


at 


( 1  +  Ky)  ax 


1 


By  (  1  +  Ky) 


Bp 


p(  1  +  Ky)  dx 
K  dv 


r  dy  r(l  +  Ky)  dx 

cos  p  dv 


a  V 


(1  +  Ky)  dxdy  r(  1  +  Ky)  dx 


K 


du  K  cos  p 


K^ 


>u  + 


(  1  +  Ky)^  dx  (  1  +  Ky)  dy  r(  1  +  Ky)  (  1  +  Ky)- 


d^u  cosp  du  1  d^u 


dr  dw 


a^w 


dy‘^  r 

(Flow  nornnal  to  body) 


dy 


dy^  r^(l  +  Ky)  dx  8y  r(i+Ky)  dxdy 


dv  dv 

_  4-  vr _  + 

u  dv  w  dv  w 

-  ^  V  I 

dt  dy 

(  1  +  Ky)  dx  r  dy  r 

1  dp 

V  \J 

'  CO  s  p  dw  1  d^  w 

+ 

p  9y 

^  K 

r-^  dy  r  dydy 

COS  p  - 


K 


1  a^v 


(1  +  Ky) 
1 


a^v 


uK 


dr  dv 


dr 


uKy 


dK 

dv 

4- 

K  du 

u 

-f. 

dK 

dx 

dx 

(1 

+  Ky)^  ax 

(  1  +  Ky) 

dK 

X  _ 

1 

a-^u 

X 

1 

dr  du 

dx 

(1 

+  Ky)  dydx 

r(l  +Ky) 

dx  dy. 

(Flow  around  body) 


dw  w  dw 

u  dw 

dw  uw 

+  V  + 

dr  vw 

CO  S  p 

dt  r  dy  (1 

+  Ky)  ax 

ay  r(l  +  Ky) 

ax  r 

1  dp 

1 

d^u  1 

dr  du 

w 

d^r 

pr  ay 

.r(l  +  Ky) 

dydx  r  ^(  1  +  Ky) 

dx  dy 

r(  1  +  Ky)^ 

dx^ 

1 

dr  dw  wy 

dK  dr 

w 

X 

dr 

r(l  +  Ky)^  ax  ax  r(  1  +  Ky)^  dx  dx  r^(  1  +  Ky)^  dx 
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1  3^w  y  dK  8w  cos  p  dw 

(1  +  Ky)^  ax^  (1  +  Ky)^  Ox  dx  r  dy 

K  cos  P  CO  8^  P  a^w  K  aw 

-  - w  +  - w  -  -  -  ■  ■ 

r(l  +  Ky)  dy^  (1  +  Ky)  dy 

1  a^v  K  av  cos  p  av 

r  ay  ay  r(l  +  Ky)  ay  ay. 

(Continuity) 

au  u  ar  av  (1  +  Ky)co8p  (1  +  Ky)  aw 

-  + - +  ( 1  +  Ky) -  +  vK  + - V  + - =  0  (20) 

ax  r  ax  ay  r  r  ay 


REYNOLDS'  EQUATIONS 


Applying  Reynolds'  procedure  (Ref.  11)  to  Eq.  20,  the  velocity  conn- 
ponents  and  pressure  for  the  turbulent  flow  are  considered  as  being 
made  up  of  a  mean  component  designated  by  a  bar  and  a  fluctuation 
component  designated  by  a  prime. 

u  =  u  +  u'  v='v  +  v'  w=w  +  w'  P  =  P  +  P'  (21) 


A  mean  component  is  determined  by  a  time-average  at  a  fixed  point  in 
space,  so  that 


"xP  =  0  =  0  w  ^  =  0  ^  =  0  (22) 

Letting  f  and  g  be  any  two  dependent  variables  and  z  be  any  inde¬ 
pendent  variable,  the  following  rules  for  operating  on  a  mean  time- 
average  apply. 


f  =  f  f  +  g=  £+  g 

_ _  ‘af  87 

fg  =  fg  —  =  — 
dz  az 


(23) 


Restricting  this  development  to  the  case  of  steady- state  flow,  and  since 
the  flow  is  axisymmetric,  the  following  relationships  hold. 

a?  aT  _ 

—  =  0  -  =  0  w  =  0  (24) 

at  ay 
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Substituting  Eq.  21  into  Eq.  20,  taking  the  mean  time-average  of  Eq.  20, 
and  applying  the  relationships  given  in  Eq.  22,  23,  and  24,  the  Reynolds 
equations  for  the  axi  symmet r i c ,  steady- state  flow  over  a  body  of  revo¬ 
lution  arc  obtained. 

(Flow  along  body) 

u  3u  0u  uv^  1  3u'^  9u’v' 

- +  V -  +  - K  + - + - 

(  1  +  Ky)  8x  ay  (  1  +  Ky)  (  1  +  Ky)  Bx  By 

ar  2u'v'  u'v'  dr 

+ - +  - K  t- -  CO  s  p - 

r(l  +  Ky)  ax  (1  +  Ky)  r  r(  1  +  Ky)  ax 

1  ap"  1  a^^  cosp  av  k  a^ 

r - V - + - - 

p(l  +  Ky)  ax  1(1  +  Ky)  axay  r(  1  +  Ky)  ax  +  Ky)'^  ax 

K  air  K  cos  p  K^  a'^rr  cos  p  air 

(1  +  Ky)  ay  r(l  +  Ky)  (1  +  Ky)^  ay^  r  dy . 
(Flow  normal  to  body) 

dv  \x  BV  av*^  1  au'v' 

V -  + - K  + -  + - 

ay  ( 1  +  Ky)  ax  (  1  +  Ky)  ay  (  1  +  Ky)  ax 

u'v'  ar  v'^  v'^  w'^ 

+ - K  + - K  + -  cosp - cosp 

r(  1  +  Ky)  ax  ( 1  +  Ky)  (  1  +  Ky)  r  r 


m.i 

1  a^  V 

1 

Br 

<1 

+ _ 

y  aK  av 

P  8y  1 

(  1  +  Ky)^  ax*^ 

r(  1  +  Ky)^ 

ax 

ax 

(1 

+  Ky)^  ax  ax 

K 

ail  u 

dK 

uK 

ar 

liKy  aK 

(1  +  Ky)^  ax  (1  +  Ky)^  dx  r(  1  +  Ky)^  ax  (1  +  Ky)^  ax 

1  a^u  1  ar  an 

+  — -f  ■  -  — 

(1  +  Ky)  ay  ax  r(l  +  Ky)  ax  By  . 
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(Flow  around  body) 


1  0u’w'  0v*w*  2u*w'  9r  v*w*  2v*w' 

- + -  + - + - K  + -  cos  p  =  0 

( 1  +  Ky)  dx  0y  r(  1  +  Ky)  0x  ( 1  +  Ky)  r 

(Continuity) 

0u  u  0r  0"^  (1  +  Ky)co8p 

-  + - +  (  1  +  Ky) -  +  vK  + - V  =  0  (25) 

0x  r  0x  0y  r 


BOUNDARY-LAYER  EQUATIONS 

In  order  to  apply  Eq,  25  in  a  practical  manner  to  the  solution  of 
physical  problems,  a  great  simplification  is  necessary.  This  can  be 
accomplished,  following  after  Prandtl  (Ref,  12),  by  considering  the 
flow  as  composed  of  two  regions:  a  region  close  to  the  body  where  the 
viscous  and  turbulent  stresses  strongly  affect  the  flow  (the  boundary 
layer)  and  a  region  outside  the  boundary  layer  where  the  flow  is  essen¬ 
tially  potential.  The  fluid  velocity  in  the  x  direction  at  the  edge  of  the 
boundary  layer  is  denoted  as  U,  and  for  elongated  bodies  of  practical 
interest  (excluding  the  region  around  the  nose  where  laminar  flow  nor¬ 
mally  exists)  there  results 

U 

-  =  0(1)  (26) 

Uo 

where  Uq  is  the  undisturbed  free- stream  velocity.  A  length  L  may  be 
defined  by  the  relationship 


For  most  practical  problems,  the  length,  L,  so  defined  will  be  of  the 
order  of  magnitude  of  the  length  of  the  body  and  will  satisfy  the  con¬ 
dition 


6 

—  <<  1  (28) 
L 

The  radius  of  curvature,  l/K,  of  the  body  surface  in  the  meridian  plane 
will  normally  be  of  the  order  of  magnitude  of  the  body  length  over  that 
portion  (excluding  the  region  around  the  nose  where  laminar  flow  nor¬ 
mally  exists)  of  practical  bodies  that  are  in  turbulent  flow.  Hence  it  is 
not  unreasonable  to  assume  that 
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(^9) 


Since  this  developnnent  must  be  applicable  to  the  case  where  the 
boundary-layer  thickness  is  not  necessarily  small  compared  to  the 
body  radiusp  it  must  be  assumed  that 

r  =  0(6)  (30) 


Applying  Eq.  26  through  30,  and  using  the  normal  order  of  magnitude 
considerations  for  determining  the  boundary  layer  equations  from  the 
Navier- Stokes  equations,  the  relative  magnitudes  of  the  terms  in  Eq,  25 
involving  the  mean  velocity  components  can  be  determined.  However, 
such  an  analysis  will  tell  nothing  about  the  order  of  magnitude  of  the 
terms  involving  the  mean  time-average  of  the  products  and  squares  of 
the  fluctuation  velocity  components.  To  obtain  some  idea  of  the  magni¬ 
tude  of  these  terms  it  is  necessary  to  resort  to  experimental  data. 

Some  of  the  most  reliable  measurements  of  the  fluctuation  velocity 
components  in  a  turbulent  boundary  layer  approaching  separation  have 
been  made  by  Newman  (Ref.  13).  Using  Eq.  27,  an  appropriate  value 
of  L  was  determined  for  the  second  series  of  tests  in  Ref,  13,  The  data 
from  this  second  series  was  then  analyzed  and  the  following  orders  of 
magnitude  determined. 
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Since  Newman’s  experiments  did  not  include  a  measurement  of  it 

was  necessary  to  examine  some  other  data.  The  experiments  of 

Ref.  14.  15,  and  16  indicate  that  the  values  of  w*^  tend  to  lie  between 

those  of  and  u’^.  Hence  it  is  reasonable  to  assume  that 


With  the  addition  of  Eq,  31  and  32,  a  complete  order  of  magnitude 
analysis  of  Eq.  25  is  possible.  In  making  such  an  analysis,  it  is  cus¬ 
tomary  to  divide  the  first  three  parts  of  Eq.  25  by  U^/L  and  the  last 
part  by  Uq/L#  putting  them  in  nondimensional  form.  Since  inertia  and 
viscous  effects  are  of  the  same  order  of  magnitude  in  the  boundary 
layer,  an  examination  of  the  inertia  and  viscous  terms  in  this  non- 
dimensional  form  shows  that 


V 


<) 


riri 


All  the  terms  will  then  be  of  the  order  of  magnitude  1,  6/L,  (6/L)^, 
(6/L)^,  (6/L)"^.  Retaining  only  those  terms  of  the  order  of  magnitude 
1,  there  results 


(Flow  along  body) 


r  a(u/Uo) 


1  iL  8(y/L) 


.(1 


1 


d 


r  u*v'\ 
L  U|  ) 


(33) 
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(Flow  normal  to  body) 


( Continuity) 


(34) 


(35) 


The  equation  for  flow  around  the  body  has  been  dropped  since  it  involves 
terms  containing  only  the  fluctuation  velocities. 

Equation  34  expresses  the  variation  of  f^^atic  pressure  across  the 
boundary  layer.  Integrating  Eq.  34  with  respect  to  y/L  and  taking  the 
partial  derivative  with  respect  to  x/L,  there  results 


Z4 
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Equation  36  may  now  be  substituted  into  Eq.  33;  however,  the  magnitude 
of  the  terms  should  be  examined  first.  The  first  term  on  the  right  side 
of  Eq.  36  expresses  the  pressure  gradient  at  the  edge  of  the  boundary 
layer  arising  from  the  external,  essentially  potential  flow.  The  re¬ 
maining  terms  express  the  variation  of  the  pressure  gradient  across 
the  boundary  layer.  An  order  of  magnitude  analysis  of  these  remaining 
terms  shows^  that  they  are  of  the  order  6/L.  Hence,  the  retaining  of 
these  remaining  terms  in  Eq.  36  when  substituting  into  Eq.  33  would  be 
inconsistent  with  the  approximation  already  made  when  terms  of  the 
order  of  b/L  were  discarded  in  arriving  at  Eq.  33.  Therefore,  within 
the  approximation  of  retaining  terms  of  the  order  of  magnitude  1  in  the 
nondimensional  Reynolds  equations,  it  may  be  assumed  that  the  pres¬ 
sure  is  constant  across  the  boundary  layer  at  the  value  given  by  the 
external  flow  at  the  edge  of  the  boundary  layer  (y  =6).  Since  the  ex¬ 
ternal  flow  is  assumed  to  be  essentially  potential,  Bernoulli's  equation 
applies,  so  that 


af_L) 

\pu|;  \pU2/y=6 


d 


u 


d 


Uo 


Substituting  the  expression  above  into  Eq.  33  and  returning  to  dimen¬ 
sional  form,  the  first-order  boundary-layer  equations  for  the  turbulent, 
axisymmetric  flow  over  a  body  of  revolution  are  obtained. 


0u  0U 

u -  +  "v - 

0x  0y 


dU  1 

=  U -  +  V  - 

d>-  r 


0y 


1  0(ru)  1  0(rv) 


- + - -  0 

r  0x  r  0y 


1  0(ru’v') 
r  0y 


(37) 


It  has  been  pointed  out  by  several  investigators  (e,  g,  ,  Ref,  13  and 
17)  that  the  above  first-order  boundary-layer  equations  in  their  inte¬ 
grated  form  (the  momentum  integral  equation)  do  not  yield  satisfactory 
results  in  a  strong  adverse  pressure  gradient  when  the  shear  stress  at 


This  result  i..  expected  since  the  terms  in  Eq.  34  are  of  the  order 
of  magnitude  1  so  that  an  integration  across  the  boundary  layer  with 
respect  to  y/L  yield  a  variation  in  pressure  across  the  boundary  of  the 
order  S/L. 
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the  wall  is  being  determined  from  a  knowledge  of  the  pressure  gradient 
and  the  velocity  profiles.  It  i  s  pointed  out  that  cer  tain  terms  of  the 
order  of  6/L  (namely,  the  turbulent  normal  stress  terms  in  the  equation 
for  flow  along  the  body  and  the  terms  expressing  the  variation  in  pres¬ 
sure  across  the  boundary  layer  in  the  equation  for  flow  normal  to  the 
body)  must  be  retained  in  the  nondimen sional  Reynolds  equation  to  yield 
anything  approaching  a  satisfactory  solution.  This  is  true,  because  in 
a  strong  adverse  pressure  gradient,  the  shear  stress  term  is  very 
small  compared  to  the  terms  expressing  the  pressure  gradient  and  the 
rate  of  change  of  the  velocity  profiles.  Hence,  the  magnitude  of  a  small 
term  in  an  equation  containing  large  terms  is  being  sought  where  a 
number  of  small  terms  have  already  been  discarded  from  the  equation. 
In  effect,  the  determination  of  a  term  of  second-order  magnitude  is 
being  attempted  using  a  first-order  equation. 

In  the  present  work,  it  is  planned  to  use  the  first-order  boundary- 
layer  equations  to  determine  the  velocity  profiles  from  a  knowledge  of 
the  pressure  gradient  and  the  shear  at  the  wall.  Thus  a  first-order 
effect  is  being  solved  for  using  a  first-order  equation,  and  there  is 
good  reason  to  believe  that  satisfactory  results  essentially  to  the  point 
of  separation  can  be  obtained. 

The  expressions  for  the  laminar  shear  stress,  Tj^,  and  the  turbulent 
shear  stress,  are 


au 

=  p_ 

ay 


(38) 


=  -pu*v* 

Letting 

T  equal  the 

total 

shear  st: 

res 3  (t 

+  T^)  and 

into  Eq. 

37,  there  is 

obtai 

ned 

au 

air 

dU 

1  1 

a(rT) 

u - 

+  V  -  = 

U -  + 

— 

ax 

ay 

dx 

P  r 

ay 

(39) 


1  a(ru)  1  a(rv) 

-  - -  + - =  0 

r  8x  r  ay 


(40) 


These  are  the  first-order  boundary- layer  equations  in  their  most  com¬ 
pact  form  for  the  axi  symmetric  flow  over  a  body  of  revolution. 
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MOMENTUM  INTEGRAL  EQUATION 


The  integration  of  Eq.  39  across  (x  =  constant)  the  boundary  layer 
(from  y  =  0  to  y  =  6)  can  be  carried  out  by  making  use  of  the  last  of 
Eq,  17  and  by  using  the  continuity  equation  40  to  eliminate  "v.  There 
re  suit  s 


where 


(41) 

(42) 

(43) 

(44) 


Equation  41  is  the  momentum  integral  equation  for  the  axi symmetric 
flow  over  a  body  of  revolution,  6^  and  0^.  must  be  looked  upon  as 
pseudo-displacement  and  pseudo- momentum  thicknesses  since  they  do 
not  retain  the  physical  significance  their  two-dimensional  counterparts 
have.  If  l/rQ  is  set  equal  to  zero  in  Eq.4l,  43,  and  44,  Eq.41  reduces 
to  the  more  familiar  two-dimensional  momentum  integral  equation  of 
vonKarman  (Ref.  18)  and  6f  and  0r  red  ce  to  the  displacement  thickness 
and  momentum  thickness  of  two-dimensional  flow. 


ENERGY  INTEGRAL  EQUATION 


Another  integral  relationship  can  be  determined  from  Eq.  39  by  first 
multiplying  Eq.  39  by  u  and  then  carrying  out  the  integration  across  the 
boundary  layer  in  a  manner  similar  to  that  used  in  obtaining  Eq,  4  1. 

The  resulting  expression  is 


d6?^  6?*  dU  6^  dro 

-  +  3 - + - 

dx  U  dx  r^  dx 

2  ^6  Sn*  2  cos  p  0u 

= - r  I  T  - dy  + - - -  I  T - ydy 

pU^  Jo  By  pU^  r^  Jo  By 


(45) 
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where 


(46) 


Equation  45  is  the  energy  integral  equation  for  the  axi  symmetric  turbu¬ 
lent  flow  over  a  body  of  revolution,  and  6r  is  a  pseudo- ene rgy  thickness 
since  it  does  not  retain  the  physical  significance  of  its  two-dimensional 
counterpart.  If  I/tq  is  set  equal  to  zero  in  Eq.  45  and  46,  Eq.  45  re¬ 
duces  to  the  two-dimensional  energy  integral  equation  first  proposed  by 
Wieghardt  (Ref.  19)  in  a  slightly  different  form  for  laminar  flow,  and 
reduces  to  the  energy  thickness  of  two-dimensional  flow.  A  con¬ 
venient  thickness  ratio  to  use  in  conjunction  with  the  energy  integral 
equation  is 

_ 

Hj.  = -  (47) 

er 


EVALUATION  OF  THE  THICKNESS  PARAMETERS 


Using  the  identities 


the  expressions  for  the  thickness  parameters,  Eq.  43,  44,  and  46  may 
be  written  in  the  following  form. 
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5 


* 

r 


V* 


6  U 


U  -  u  /y\  V*  6  /•!  U  -  u  y  /yl 

- d-  + - cosp  I - d- 

\6/  U  Jo  6  \6/ 


U-u/y\  v^6  /•lU-uy/yl 

- dl—  )  + - cosp  I - “  1 

16/  U  Tq  Jo  6  \6/ 


/•I  /  U  -  u  /y  \  /  6  /•' 

I  (  -  I  —  1  ”  I  -  J  -  COS  P  I 

•/o  \  /  \  6  /  \  u  /  Tq  -/q 


6 


■X-* 

r 


6 


+  2 


6 


cos  P 


Using  the  analytical  approximation,  Eq,  14,  to  Coles'  two-parameter 
mean  velocity  profiles,  the  integrals  in  the  expressions  above  may  be 
evaluated  (see  Appendix).  The  expressions  for  Hj.,  Hj.,  and  6  then 
become 
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Hr 


6 

N  j  +  AN2  +  cos 

_ 

6  Vn  1  f 

Ni  +  ANi  + -  cosp(N3  +  AN4) - N5  +  AN 

To  U  «  I 


H 


r 


16  1  *  ?  }  h 

Ni  +  AN2+ cosP(N3+AN4) - N5  +  AN^  +  A^N7  + — cos  P(N8  ^  AN9  •I' A^N  iq)  ^  - 

^o  J  U  «  I  Fq  J  \  ' 


6  V*  1  ’ 

N  1  +  ANa  + - cos  ^(N3  +  N4) - Ns  +  AN^ 

r©  U  «  I 


Or 

6  - - 

V,  1  r  6  1  /  *  f 

- Ni  +  AN2  + -  co»P(N3  +  N4)  -I - |  - INi 

U  «  I  To  J  \  U  /  «^  I 


where  the  numerical  constants,  Nj  through  Nig,  are  given  in  the  Appe 
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.•P(N3  +  AN4) 

6 

Hf,  +  a2n7  + —  coap(N8  +  AN9  +  A^Nio) 
'o 


(48) 


Ni  1  +  AN12  +  A^Ni3+  A^Ni4^— cos  pCNjj  +  AN  ifc  +  A^N  17  +  A^N  jg) 
6 

6  +  A^N7  +  —  cos  P(N8'«-  AN9  +  A^N  iq) 

>■0 


N5  +  ANfc  +  A^Nf  + - (Ng  +  AN9  +  A^Nio) 


>endix. 
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THE  SHEAR  STRESS  AT  THE  WALL 

The  shear  stress  at  the  wall,  Tq  =  pv^  ,i8  specified  if  the  velocity 
at  the  edge  of  the  boundary,  U,  the  boundary* layer  thickness,  6,  and 
the  profile  parameter.  A,  are  known.  This  is  seen  by  applying  the 
boundary  condition'll  =  U  when  y  ==  6  to  Eq*  13. 

U  1  6v^  If  1  /  n  \1 

—  =  «  In  -  +  C- - +  Ajl - I 

K  V  iclm  \  m/J 

This  equation  can  be  put  in  a  more  useful  form  in  which  the  boundary 
layer  thickness,  6,  plays  a  secondary  role  by  using  Eq.  50.  The  follow¬ 
ing  implicit  equation  for  /U  results. 

U  1  OrU  1 

-  =  —  In  -  +  C  +  — 

ic  V  ,c 


1 

”  \ 

- +  A  1 

1 

,  m 

k  ml, 

1  1 
-  In  - 
K  Is 


U 


Ni  +  ANz  + -  cos  P(N3  +  AN4) 


N5  +  AN^  +  A^N7  + -  cos  P(Na  +  AN9  +  A^ 


'Nio)]} 


(51) 


DETERMINATION  OF  THE  ENERGY  DISSIPATION 
The  two  integrals 


T 


9u  Z  cos  p 

- dy  + - -  - 

dy  pU^  Tq 


8u  2  1  ^  an 

T  - ydy  - - -  -  I  — dy 

8y  pU-*  Tq  Jo  By 


(52) 


on  the  right  side  of  the  energy  integral  equation  (45),  represent  the  rate 
at  which  organized  kinetic  energy  in  the  boundary  layer  is  dissipated  in¬ 
to  random  kinetic  energy  (turbulence)  and  heat.  In  order  to  evaluate 
these  integrals,  a  knowledge  of  the  variation  of  the  shear  stress,  r, 
across  the  boundary  layer  is  needed.  This  variation  can  be  found  by 
using  a  momentum  integral  equation  in  which  the  integration  is  not 
carried  completely  across  the  boundary  layer  (y  =  0  to  y  =  6)  but  rather 
partially  across  the  boundary  layer  (y  =  0  to  y  =  y).  Using  the  first- 
order  boundary-layer  equations,  39  and  40,  the  incompletely  integrated 
momentum  integral  equation  takes  the  form 
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1  1  _ 

-  rr  =  -  ToTo  ’  “  /  -  dy 

p  p  ./o  ^x 


J^y  d{  ru)  0u  du  /-y 

^  u  -  dy  +1  ru  —  dy  -  U  —  I  rdy  (53) 

o  Jo  ^  dx  «/o 

This  result  reduces  to  the  momentum  integral  equation  41,  if  evaluated 
at  y  -  6. 

The  question  arises  at  this  point  whether  the  first-order  equation  53 
can  yield  a  reasonably  good  solution  for  the  variation  of  t  across  the 
boundary  layer  or  whether  second-order  effects  must  be  included. 

Coles  (Ref.  1)  compared  the  variation  in  t  across  the  boundary  layer 
given  by  the  two-dimensional  form®  of  Eq.  53  with  Newman's  second 
series  of  tests  (Ref.  13).  The  integrals  in  the  two-dimensional  form  of 
Eq.  53  were  evaluated  using  the  mean  velocity  profiles  given  by  Eq.  5 
which  had  been  fit  to  the  experimental  mean  velocity  profiles.  Coles 
found  that  the  general  behavior  of  t  versus  y  was  reasonably  well  pre¬ 
dicted  but  that  both  the  boundary  condition  t  =  when  y  =  0  and  t  =  0 
when  y  =  6  could  not  always  be  satisfied  since  the  experimental  data 
did  not  conform  to  the  momentum  integral  equation  at  all  points  in  the 
flow.  This  noncomformity  was  attributed  primarily  to  a  deviation  of 
the  flow  from  true  two-dimensionality.  By  a  small  alteration  in  the 
continuity  equation.  Coles  made  the  calculated  values  of  t  satisfy  both 
boundary  conditions  and  brought  them  in  good  agreement  with  the  ex¬ 
perimental  values.  Since  it  is  planned  to  solve  the  momentum  integral 
equation  simultaneously  with  the  energy  integral  equation,  Eq.  53  will 
satisfy  both  the  boundary  conditions  t  =  Tq  when  y  =  0  and  t  =  0  when 
y  =  and  based  on  Coles'  findings  should  give  a  satisfactory  solution 
for  T  versus  y. 

Substituting  Eq.  53  into  Eq.  52,  the  energy  dissipation  integrals  may 
be  written,  dropping  the  multiplier  2/U^, 


1  BU. 

P  Jo  By 


1  cos  p 
dy  +  —  —  -  — 


8y  p 

ro 

1  /-< 

i  8u 

=  -  >^0^0  I 

p  Jo 

8y 

^6  /  8u 

+  1  —  i 

f  Ti 

Xo  du 

dy 
5  i 


ydy 

'  du  /•y  d(ru) 


dx 


•) 


dy  Idy 


8<ru) 


8y 


I  0u  _  8u  \ 

dy  +  I  I —  I  ru - dyldy 

Jo  \  dy  Jo  8x  / 

dU  I  8u  /•y  \ 

-U -  I  I -  f  rdyjdy  (54) 

dx  Jo  \  dy  Jo  I 

The  two-dimensional  form  of  Eq.  53  is  obtained  by  dividing  through 
by  r  and  letting  ro 
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Using  the  last  of  Fq.  17  and  integration  by  parts,  Eq.  54  can  be  expanded 
into  a  form  containing  much  simpler  integrals.  In  this  operation,  one 
approximation  is  indicated  based  on  the  previous  order  of  magnitude 
analysis.  Differentiating  the  last  of  Eq.  17  with  respect  to  x,  there 
results 

dr  dro  dp 

—  = -  -  y  sin  p - 

dx  dx  dx 

But  from  Fig.  3  it  is  seen  that 

dr^  dp 

sin  p  =  ■  —  -  =  -K 

dx  dx 

Therefore, 

dr  dro 

—  = - (1  +  yK) 

dx  dx 

But  ymax  ”  ^  0(  l/L»)  so  that  yK  =  0(6/Li),  and  since  6/L  <<  1 

it  follows  that 

ar  dro 

—  ar - 

dx  dx 

Equation  54  then  takes  on  the  following  form  in  which  the  multiplier 
has  been  reintroduced  and  y  has  been  nondimensionali zed  by  6. 
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6  6  /•lu8uy/y\  6  dU  /•!  /  y  \ 

4_cosP-  / - d-  -2 - /  d- 

rj,  UJbU8x6\6/  UdxJo  I 

6  du  6  /*^  y  /  y  \  ®  “  /  y  \ 

-2 - cosp  I  _d  -l+2 - /  — d  - 

U  dx  ro  Jo  6  \s  I  Udx  Jo  u  \6/ 

6  dU  6  "S  y  /  y  \ 

+  2 - cosp  I - ^(~  I  (55) 

U  dx  To  Jo  U  6  \6/ 


The  8Cr/8x  is  readily  determined  from  Eq.  13 
8u  u  U  dv. 


8x  U  Vo  dx 


1  dv^  IdA  An  dSWyV” 

x  dx  K  \  dx  6  dx  /  \  6  / 

IKI 


n  dA  (1  +  An)  d6" 
m  dx  6  dxj 


Using  £q.  56  and  the  following  identities. 


u 

U 


v^/U  -  u 


U 


V* 


-1  =1-2 
U; 


er--T(=^")-e)ra‘-s‘(=;r7 


(56) 


”  '  V*  /  \  U  /  \  V,  /  \  U  /  \  V* 

all  the  integrals  on  the  right  side  of  Eq.  55  may  be  put  into  the  form 


ri^'rer-ei 


(q  =  integer) 


These  integrals  may  be  evaluated  using  Eq.  14  (see  Appendix)  and  yield 
the  final  expression  for  the  energy  dissipation  integrals. 
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n 


m 


m 


+ -  cospl  —  N36  -  N30  +  A  — N37  -  N31)  +  A^l  —  N38  -  N32 


m 


m 


d6  1 

+  2 - 

dx  U  « 


—  11 
U  xl 


N23  +  A(N24  +  nN23  -  nNi9)  +  A^(nN24  -  nN2o) 


+  —  co8p[N25  +  A(N26  +  *'^25  -  nN2i)  +"^^(^^26  -  nN22)l 
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3/v»v2  1  I 

+  — I  — I  — |-N33  +  A(nN27  “  "N33  -  N34) 


+  A2(nN28  -  nN34  -  N35)  +  A3(nN29  -  nN35) 


+ -  cosp[-N36  +  A(nN3o  "  nN36  -  N37) 

*•0 

+  A2(nN3|  -  nN3Y  -  N30)  +  A^(nN32  -  nN3g)] 

6  d^o  V#  1 


Tq  dx  U  « 


V4  1 

- (_N5  -  AN6  -  A^N7) 

U  K 


+ 


2  1 

I  —  (Nil  +  AN12  -!  a2Ni3  +  a3Ni4) 


+ 


5  dU  V*  1 
U  dx  U  « 


6 

.Ni  -  AN2  + —  cosp(-N3  -  AN4) 


(57) 


where  the  numerical  constants  Nj  through  N38  are  given  in  the  Appendix 
and  n  and  m  are  regarded  as  equal  to  2.50  and  2.75,  as  mentioned  pre¬ 
viously.  Also  the  substitution  Tq/pU^  =  (v^/U)^  has  been  made. 


THE  BOUNDARY-LAYER  SOLUTION 

Equations  41,  42,  45,  47,  48,  49 »  50,  51,  and  57  constitute  a  set  of 
nine  differential,  algebraic,  and  transcendental  equations  in  the  nine 
unknowns: 
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It  i  s  assumed  the  pressure  distribution  at  the  outer  edge  of  the  boundary 
layer  U  =  f(x)  and  the  body  shape  l/r^  =  f(x)  are  known.  A  solution  of 
these  equations  for  l/rQ  /  0  is  a  first-order  solution  lor  the  turbulent 
boundary  layer  on  a  body  of  revolution  in  axisvmmetric  flow.  For 
l/r^  =  0,  it  is  a  first-order  solution  for  the  turbulent  flow  on  a  two- 
dimensional  body.  As  brought  out  earlier  in  the  development,  both 
solutions  are  restricted  to  the  case  where  the  radius  of  curvature,  l/K, 
of  the  body  surface  in  the  direction  of  flow  is  large  compared  to  the 
boundary-layer  thickness. 

An  iterative  method  for  obtaining  a  simultaneous  solution  of  these 
nine  equations  is  under  development,  and  programming  of  the  iterative 
method  on  the  IBM  7090  will  be  undertaken  in  the  near  future. 

To  commence  such  an  iterative  solution,  starting  values  of  ,  H^., 
and  X  must  be  known.  Starting  values  of  0r  and  x  may  be  determined 
from  a  laminar  boundary-layer  solution  from  the  stagnation  point  to 
the  transition  point.  The  Pohlhausen  method  (Ref.  20)  offers  a  conven¬ 
ient  way  of  obtaining  the  laminar  solution  since  it  may  be  applied  in  a 
simple  approximate  form  to  both  the  two-dimensional  body  (Ref.  2  1) 
and  the  body  of  revolution  (Ref.  22).  A  reliable  starting  value  of  Hj.  is 
not  easily  obtained  because  of  the  scarcity  of  information  on  the  change 
of  this  thickness  ratio  through  transition,  particularly  in  a  pressure 
gradient.  It  is  necessary,  therefore,  to  determine  a  starting  value  of 
Hr  by  some  rather  uncertain  means  such  as  is  suggested  in  Ref,  23  and 
24.  This  is  admittedly  a  weak  link  in  the  boundary- layer  solution. 


THE  DISPLACEMENT  THICKNESS  ON  A  BODY 
OF  REVOLUTION 

For  the  axisymmetric  flow  on  a  body  of  revolution,  6r  is  a  pseudo¬ 
displacement  thickne  ss ,  as  mentioned  earlier.  It  does  not  represent 
the  increase  in  the  body  dimension  (normal  to  the  surfaceY  that  would 
displace  the  streamlines  of  a  totally  potential  flow  around  the  body  by 
the  same  amount  they  are  displaced  as  a  result  of  the  retardation  of 
the  flow  in  the  boundary  layer.  The  true  displacement  thickness,  6^, 
which  does  represent  this  condition,  is  determined,  as  pointed  out  by 
Granville  (Ref.  25)  by  the  relation 


6* 


co  s  (36r 


1 

—  co  S  P 


In  the  limit  as  l/r^  0,  the  equation  above  yields  6*  =  6^,  which  is  the 
expected  result  from  the  definition  of  6^* 
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It  is  often  customary  to  compute  the  potential  pressure  distribution 
around  a  body  increased  in  dimension  by  the  displacement  thickness. 
This  is  done  to  account  for,  in  part,  the  effect  of  the  boundary  layer  on 
the  potential  flow  outside  the  boundary  layer.  To  carry  out  this  calcu¬ 
lation,  the  values  of  6'^  are  needed. 


DISCUSSION  AND  CONCLUSIONS 

Although  most  of  the  development  leading  from  the  Navier-Stokes 
equations  to  the  momentum  and  energy  integral  equations  can  be  found 
in  various  references,  it  was  included  in  this  report  in  order  to  empha¬ 
size  the  assumptions  and  restrictions  underlying  the  first-order 
boundary-layer  equations  and  to  give  a  line  of  unbroken  continuity  to 
the  development. 

The  foregoing  development  has  shown  that  the  use  of  Coles*  univer¬ 
sal,  two-parameter,  mean  velocity  profiles  made  possible  a  scheme 
for  obtaining  an  essentially  exact  solution  of  the  first-order  turbulent 
boundary- layer  equations  in  integrated  form.  This  was  possible  be¬ 
cause  these  profiles  inherently  contain  sufficient  information  to  deter¬ 
mine  not  only  the  shear  stress  at  the  wall  but  also  the  variation  of 
shea'  stress  across  the  boundary  layer,  making  possible  the  determi¬ 
nation  of  the  rate  of  energy  dissipation  in  the  boundary  layer.  Due  to 
tne  fairly  rigorous  development  of  this  method  and  the  fact  that  velocity 
profiles  can  be  accurately  represented  essentially  to  the  point  of  sepa¬ 
ration,  it  is  anticipated  that  an  improved  technique  for  estimating  the 
inflow  velocity  profiles  to  propellers  and  the  skin-friction  drag  of 
bodies  of  revolution  will  result. 
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Appendix 

EVALUATION  OF  INTEGRALS 

The  development  in  the  body  of  this  report  reqiires  the  evaluation 
of  integrals  of  the  form 


using  Eq.  14  (repeated  below) 

U  -  u  1  y  A 

-  -r  .  —  In  —  +  —  I 

#c  6  K  , 


(q  integ<>r) 


(  H) 


where  n  ~  and  m  =  2.75.  These  intt'grations  are  straiglitforwa  rd 

but  they  are  rather  long  and  tedious  when  q  is  greater  than  one*.  In 
carrying  out  these  integrations,  the  terms  involving  ln(y/!>)  lead  to 
improper  integrals  since  Eq.  14  does  not  properly  dt‘scribe  the  flow  in 
the  region  very  near  the  wall.  These  improper  inti*grals  are,  how¬ 
ever,  convergent,  and  the  effect  on  the  value  of  these  intt'grals  due  to 
tile  deviation  from  the  true  flow  in  the  transition  zone  and  laminar  sub¬ 
layer  is  not  significant. 


K 


-  N  1  4-  AN^ 


K 


-  N3  +  AN4 


N5  +  AN^  + 


-  Ng  +  ANq  +  A‘“N  jq 
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^L{ 

^L'i 

I 


1  /  U  -  u 


d(  -  I  =  Ni  1  +  AN  12  +  A^Ni3  +  A^Ni4 


1  / U  -  u 


y  /y ' 


-d  -  =  Njs  +  ANi6  +  A^Ni7  +  A^N 

6  \b 


1  U  -  u  /  y  \2.50 


i 


o  V  ^  \  6 

1  U  -  ■u/y\3.50  I  y 


d  -  I  =  N  19  +  AN^O 
6 


£ 


’o  V  ^  \  6 

1  U  -  u/y  \2.75 


d|-  I  =  N21  +  AN22 
5 


o 


Vu  \  6 


d  -  =  N23  +  AN24 

,6 


L 


U  -  u/y\3.75  /y 

«  I  - 1-1  d  -  I  =  N25  +  AN^e 


o  V* 


Vw.  \  6 


£{ 


1  /  U  -  u\2/  y  i2.50 


dj-  I  =  N27  +  AN28  +  a2N29 

i  6 


N30  +  AN31  +  a2N32 


N33  +  AN34  +  a2n35 


N38  +  AN37  A^N38 


where 


18 
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Ni  =  0.733333334 
N2  =  0.047619048 
N3  =  0.144736842 
N4  =  0.014619884 
N5  ^  1.406495726 
=  0.116197456 
N7  ^  0.003296703 
Ng  ^  0.  U8485687 
N9  =  0.0  18400596 
N|o  =  0.000793209 
Nil  ^  4.182609023 
N|2  =  0.36895053 
N23  =  0.01447839 
N  24  ^  0,00025215 
N  25  =  0.  185477054 
N  iQ  -■■■  0.0302267  15 
N'27  =  0.001986392 
N28  =  0.0000i?1483 
N29  -  0.035918367 

These  values  are  doubtful  in  the  1 


N20  =  0.004761904 
N21  =  0.0  18731375 
N22  =  0.002736728 
N23  =  0.030085470 
N24  =  0.004102564 
N25  =  0.0  1625  1  154 
N26  ^  0.002419842 
N27  =  0,0  16054940 
N28  =  0.003292518 
N29  ^  0.000186741 
N30  =  0.005999780 
N32  =  0.001405893 
N32  =  0.000088639 
N33  =  0.0122927  19 
N34  =  0.002619501 
N35  =  0.000153119 
N36  ■=  0.004834318 
N37  =  0.00  1  163  109 
N38  =  0.0000749  1 1 

St  one  or  two  significant  digits. 
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IN  REPLY  REFER  TO: 

751/MK:ss 
4  June  1965 


From:  Commander ,  U.  S.  Naval  Ordnance  Test  Station 

To;  Distribution  List  of  NAVWEPS  Report  8510 


SubJ ;  NAVWEPS  Report  8510,  A  Turbulent  Bound ary -Layer  Calculation 
Method  Based  on  the  Law  of  the  Wall  and  the  Law  of  the  Wake, 
dated  November  1964 ;  transmittal  of  addendum  for 

Enel:  (1)  Addendum  dated  5  May  1965 


1*  It  Is  requested  that  the  discussion  provided  on  the  attached 
enclosure  be  Incorporated  In  NAVWEPS  Report  8510 « 


C.  E.  VAN  HAI^ 
By  dlrectloij^/ 
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NAVWEPS  REPORT  8510 


ADDENDUM 


It  should  be  pointed  out  that  an  additional  condition  must  be  imposed 
on  the  turbulent  boundary- layer  solution.  Experimental  evidence  shows 
that  for  flow  on  a  flat  plate  with  zero  pressure  gradient,  the  profile 
parameter  IT  in  Eq.  5  or  the  profile  parameter  A  in  Eq.  13  remains 
essentially  a  constant.  This  fact  was  first  noted  by  Schultz- Grunow 

dA 

(Ref,  3).  Therefore,  the  solution  must  yield  =  0  for  A  =  Flat  plate 
value,  =  0,  and  =  0.  (The  flow  on  a  cylinder  of  constant  radius 

is  similar  to  that  on  a  flat  plate  if  «  1.  )  This  requirement  of  the 

solution  can  be  obtained  by  adding  a  correction  term  to  the  energy  dis¬ 
sipation  function. 

A  discussion  of  this  point  had  been  postponed  for  a  future  report 
since  it  was  not  certain  in  what  form  the  correction  term  would  be  in¬ 
corporated  into  the  solution,  and  the  computer  program  that  would  be 
required  to  investigate  the  various  possibilities  was  not  available  when 
NAVWEPS  Report  8510  was  published.  However,  it  was  completely 
overlooked  that  the  application  of  this  requirement  to  the  solution  was 
not  only  desirable  but  was  essential  to  make  the  equations  independent. 
The  fact  that  the  equations  mentioned  on  page  36  contain  a  redundancy 
was  pointed  out  in  a  communication  from  Dr.  G.  E.  Gadd  of  the  National 
Physical  Laboratory,  Teddington,  England.  Since  the  variation  of  t 
across  the  boundary  layer  used  in  the  energy  integral  equation  was 
determined  from  the  partially  integrated  momentum  equation,  the 
momentum  and  energy  integral  equations  are  not  independent.  A  more 
complete  discussion  of  this  area  will  appear  in  a  forthcoming  report 
that  will  present  the  method  used  to  solve  the  equations  and  comparisons 
of  computer  solutions  with  experimental  data. 


5  May  1965 
Enclosure  (1) 


